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Abstract 

The most salient features observed around a metamagnetic transition in Sr3Ru207 are well 
captured in a simple model for spontaneous Fermi surface symmetry breaking under a magnetic 
field, without invoking a putative quantum critical point. The Fermi surface symmetry breaking 
happens in both a majority and a minority spin band but with a different magnitude of the 
order parameter, when either band is tuned close to van Hove filling by the magnetic field. The 
transition is second order for high temperature T and changes into first order for low T. The first 
order transition is accompanied by a metamagnetic transition. The uniform magnetic susceptibility 
and the specific heat coefficient show strong T dependence, especially a logT divergence at van 
Hove filling. The Fermi surface instability then cuts off such non- Fermi liquid behavior and gives 
rise to a cusp in the susceptibility and a specific heat jump at the transition temperature. 

PACS numbers: 71.18.+y, 75.30.Kz, 74.70.Pq, 71.10.Fd 
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Usually the symmetry of the Fermi surface (FS) satisfies the point-group symmetry of 
the underlying lattice structure. However, recently a symmetry-breaking Fermi surface 
deformation with a <i-wave order parameter, where the FS expands along the k x direction 
and shrinks along the k„ direction, or vice versa, was discussed in the two typical models for 

fin nn 

cuprate superconductors, the t-J[l|, [2] and Hubbard model|3|, |4j] on a square lattice. This d- 
wave type Fermi surface deformation (<iFSD) is often called <i-wave Pomeranchuk instability, 
referring to Pomeranchuk's stability criterion for isotropic Fermi liquids. However, the 
<fFSD often takes place even without violating such a criterion, since the <iFSD instability 
is usually first order for low temperature, jf], Q] and it can also happen even for strongly 
correlated systems such as those described by the t-J model. The g?FSD instability 
is driven by forward scattering processes of electrons close to the van Hove points in the 
two-dimensional Brillouin zone. As a result, the symmetry of the electronic state is reduced 
from Cfa to Cm while the lattice still retains C± v symmetry as long as no electron-phonon 
coupling is considered. The c£FSD state has the same symmetry as the so-called electronic 
nematic order, which was widely discussed in the context of charge stripes for cuprates.Q] 

Recently such a novel FS instability was reported for the ruthenate compound Sr3Ru207 
in the presence of an external magnetic field. |9( Sr 3 Ru20y is a bilayered material with metallic 
RuC>2 planes where the Ru ions form a square lattice. Ah initio calculations 10| showed that 
the electronic band structure is similar to that for the single layered ruthenate Sr 2 Ru04, 
a well-known spin-triplet superconductor. [11] The ground state of Sr 3 Ru207 is, however, 
paramagnetic. 12| By applying a magnetic field h, Sr 3 Ru207 shows a metamagnetic tran- 
sition at h — h c , around which non- Fermi liquid behavior is observed in various quanti- 
ties: resistivity, Q] specific heat,Q, Q, 16] and thermal expansion. 1?}] This non- Fermi 



liquid behavior was discussed in terms of a putative metamagnetic quantum critical end 
point (QCEP) and hence Sr 3 Ru20y was often referred to as a system with a metamag- 
netic QCEP. 18, Q, 2(3] However, subsequent experiments |s| for ultrapure crystals showed 
that the hypothetical QCEP was hidden by a dome-shaped phase transition line; a first 
order transition was confirmed at the edges of the transition line and was accompanied by a 
metamagnetic transition, while a second order transition was inferred for high temperature 
T. Grigera et a/.j^] associated this instability with a spontaneous c£FSD. Later Kee and 



Kim 



21| demonstrated a metamagnetic transition due to a first-order (iFSD phase transi- 



tion in a phenomenological model, which they discussed might be relevant to Sr 3 Ru207. A 
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direct experimental evidence of the <fFSD, however, has not been obtained so far. Other see 



narios, namely, microscopic phase separation due to Coulomb interaction [22j and magnetic 
domain formation due to long-range dipolar interactions, [23| were also proposed to explain 
the experimental data for S^R^C^. 

In this Letter, in addition to confirming a metamagnetic transition due to a first-order 



c/FSD phase transition as reported by Kee and Kim,[2l| we show that the experimental phase 
diagram for Sr 3 Ru207, and the T dependences of the uniform magnetic susceptibility and 
the specific heat are well captured in terms of the <iFSD instability near the van Hove filling, 
without invoking a putative QCEP. We analyze a simple model with a pure forward scat- 
tering interaction driving the c£FSD instability in the presence of a magnetic field. We find 
that when either the majority or minority band is tuned to van Hove filling by the magnetic 
field, the c?FSD instability occurs in both bands, but with a different magnitude of the order 
parameter. The transition is second order for high T and changes into first order for low T. 
The first order transition is accompanied by a metamagnetic transition. Both the magnetic 
susceptibility and the specific heat coefficient show strong T dependences, especially a log T 
divergence at van Hove filling. This non- Fermi liquid behavior originates from the van Hove 
singularity in the density of states of the bare dispersion. The <iFSD instability then cuts 
off such non-Fermi liquid behavior and produces a cusp in the susceptibility and a specific 
heat jump at the transition temperature. 

We investigate the ciFSD instability in the presence of a magnetic field on a square lattice. 
The minimal model reads 

h = y,&-v)<+-^ E /i* (i) 

k,cr k,o-,k,o-' k,cr 

where n£ = c k(T Ck CT counts the electron number with momentum k and spin a; c k(7 (ck CT ) is 
an electron creation (annihilation) operator; [i is the chemical potential; N is the number 
of lattice sites; h is the magnetic field. For hopping amplitudes t and t' between nearest 
and next-nearest neighbors on the square lattice, respectively, the bare dispersion relation 
is given by 

e k = — 2t(cos k x + cos k y ) — At' cos k x cos k y . (2) 
The forward scattering interaction driving the spontaneous <iFSD has the form 

iW = -g <4<4' , (3) 



with a coupling constant g > and a <i-wave form factor d^ = cosk x — cosk y . This 
ansatz mimics the structure of the effective interaction in the forward scattering channel as 
obtained for the and Hubbard model. Q| The model (1) without the magnetic field 



was extensively studied in Refs. 



and H. 

We decouple the interaction by introducing a spin-dependent mean field vf = 
~n S^k ^k(^k) an d obtain a renormalized band dispersion = e^+r/dk— \i c with 77 = J2 a V a i 
the mean fields are determined by minimizing the free energy, which is a even function with 
respect to r], and a solution of r\ > is considered. The a-summed mean filed r\ enters ££, 
and thus a finite if in general induces a finite f]~ a . The Zeeman field is absorbed completely 
in the effective chemical potential \x a = ll + ah. Since our Hamiltonian (pQ) does not allow 
momentum transfer, the mean-field theory solves our model exactly in the thermodynamic 
limit. 

The magnetization is given by m = a a(n^) = ■h a c/(Ck); an d thus the uniform 

magnetic susceptibility x is 

dm 1 / . . 

k,cr 

where f(x) = l/(e x ^ T + 1) is the Fermi function and /' is its first derivative. The electronic 
specific heat coefficient 7 = C/T can be obtained straightforwardly by the second derivative 
of the free energy with respect to T at fixed /a, 
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Tjvr y^(£k) 2 /'(£k 



T 2 N 

k,cr 



The second term Eq. (jSJ) is zero above T c and leads to a specific heat jump at T c . 

Band structure calculations for Sr 3 Ru 2 07[lo| yield 6 Fermi surfaces for h — 0. We focus 
on a FS closest to k = (ir, 0) and (0, 7r) since the (iFSD instability is driven by electrons 
near the van Hove points on a square lattice; we mimic such a FS by choosing t'/t = 0.35. 
We fix Li/t — 1 and take g/t = 1 for numerical convenience, although Sr 3 Ru207 is expected 
to have a much smaller g as we discuss later. We choose t as a unit of energy so that t = 1 
in this paper. Since the results are symmetric with respect to h — > — h and cr — » — a, we 
consider only the case h > 0. 

Figure [1] shows a phase diagram in the plane of applied magnetic field h and temperature 



T. The <iFSD transition occurs around the van Hove energy of the up-spin band {h = 0.4) 
with a second order transition for high T and a first order one for low T; end points of the 
second order line are tricritical points. Figure [2(a) shows the h dependence of the order 
parameter rj, together with r] a , at low T. Both ry and rj^ show a jump at the first order 
transition point, but with a different magnitude. The magnetization m also shows a jump 
at the first order phase transition [Fig. [2(b)]. Its upward jump with increasing h is due 
to a generic consequence of the concavity of the grand canonical potential as a function of 
h. Hence the first order transition of the rfFSD instability is necessarily accompanied by a 
metamagnetic transition. The FSs at low T are shown in Figs. [2(c) and (d) for h = 0.3 and 
0.5, respectively. The gray lines are FSs for g = and an outer (inner) FS corresponds to 
the up-spin (down-spin) band; the splitting of these FSs is due to the Zeeman energy. The 
FS instability drives a deformation of both FSs and typically leads to an open outer FS. 

Figure [3] shows the T dependence of x [Eq. Q] for several choices of h for g = 0. For 
a small field, x has a weak T dependence with a broad maximum at relatively high T and 
becomes constant for low T, i.e. Pauli paramagnetic behavior. As h moves close to the van 
Hove energy (h = 0.4), x starts to have a strong T dependence and forms a pronounced 
peak at low temperature. The peak position is pushed down to zero temperature at the van 
Hove energy, where a logT divergence appears. Similar behavior is seen when the magnetic 
field is reduced from a large h to the van Hove energy [inset of Fig. [3](a)]. Defining T* as 
the peak position of Xi we thus obtain the V-shaped T* line shown in Fig. [1] in particular 
T* goes to zero at the van Hove energy. This behavior is due to the van Hove singularity 
of the up-spin band, not due to an underlying quantum critical point. Around van Hove 
filling, the dFSD instability occurs and produces a cusp in the T dependence of x as shown 
in Fig. [3^b). A T* line is thus not defined inside the symmetry broken phase and the thin 
dashed line in Fig. [1] represents T* in the absence of the <iFSD. 

Figure [T] is very similar to the phase diagram reported by Grigera et al. jol with a first 
order transition for low T, a second order transition for high T, and a V-shaped T* line, 
which crosses the transition line near the tricritical points. In experiments, the T* line was 



determined by thermal expansion measurements. [2J] While peak positions of the magnetic 
susceptibility were not studied systematically in experiments, 24| we expect that the mag- 
netic susceptibility shows a similar V-shaped T* line around the van Hove energy. On the 
other hand, a strong T dependence of the magnetic susceptibility [Fig. Eta)] was observed in 
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experiments; 25|] the cusp of x shown in Fig. [3](b) can be tested. Experimental energy scales 
are, however, much smaller than in our result. To obtain a comparable T c ~ IK, the cou- 
pling constant g should be reduced substantially. As clarified in Ref. in the weak coupling 
limit the dome-shaped transition line of the (iFSD is characterized by a single energy scale 
CAe" 1 /*- 2 ^, where e\ is a cutoff energy and ~g = 2m*g/ir 2 is a dimensionless coupling constant 
with the effective mass m* near the van Hove energy. The result in the weak coupling limit 
was checked to be applicable even for finite g with good accuracy. [3] Therefore the phase 
diagram in Fig. [T]does not change essentially for a smaller g, which just reduces the energy 
scale of the <iFSD transition line around the van Hove point. The relative position of T* 
and T c tri also does not change appreciably, since T* has linear dependence of h around the 
van Hove energy. 

The electronic specific heat is calculated exactly in our model [see Eq. fl5])]. We first show 
7 for g = in Fig. H] as a function of T for several choices of h. For a relatively small field, 7 
shows a peak structure for high T (in a logarithmic scale), but becomes constant for lower T 
as expected for the normal Fermi liquid. As h moves close to the van Hove energy (h = 0.4), 
the peak position is shifted to lower T and a log T divergence appears at the van Hove energy. 
A similar behavior is also seen when h is reduced from a large value to the van Hove energy 
[inset in Fig. Ufa)]. For a finite g, the <iFSD instability produces a jump in the specific heat 
at the transition temperature as shown in Fig. HKb). This jump typically becomes larger as 
h is closer to the van Hove energy. Unlike the situation in the BCS theory, the ratio of the 
magnitude of the jump A7 and the normal state specific heat 7„ at T c is not a universal 
value. In particular, in the weak coupling limit T c scales as e^e -1 '^' as discussed above. 
Since 7 shows a logT dependence at van Hove filling, we obtain 7„ ~ logT c ~ g^ 1 oc g^ 1 , 
while A7 is a certain finite value. 

A T dependence of 7 very similar to Fig. Ufa) was actually obtained in experiments. [14I . 
[l^Q] The specific heat jump is not observed in experiments, but may require a more precise 
measurement, which would provide a definite evidence of the second order transition; the 

u 

authors in Ref. |9| inferred a second order transition from the magnetization and the resistivity 
measurements. 

We have shown that the most salient features observed in S^R^C^ are well captured 
in terms of the c/FSD instability near van Hove filling. In particular, the non-Fermi liquid 
behavior reported for the T dependence of \ and 7 can be associated with the van Hove 



6 



it is known 



singularity, not with a putative QCEP as usually discussed. Q Around van Hove filling, 



that various other ordering tendencies develop and probably compete 



with the <iFSD tendency. We may thus allow other interactions such as ferromagnetism and 
antiferromagnetism in our model to explore more detailed comparison with experimental 
data as well as the interplay among various ordering tendencies in the presence of the 
magnetic field. 



resistivity 
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?he most interesting future issue is how the anomalous T dependence of the 
14| observed in Sr3Ru207 above T c can be understood. A crucial question is 
whether classical <iFSD fluctuations and the van Hove singularity are sufficient to capture 
the resistivity data or whether quantum critical fluctuations originating from some QCEP 
are necessary. In the latter case, how is the QCEP related to the <iFSD instability and 
the van Hove singularity that we have discussed in the present paper? Both the chemical 
potential and the magnetic field have to be fine-tuned to reach a QCEP while it is sufficient 
for either of them to be tuned to realize van Hove filling. In this sense, a QCEP scenario 
imposes an additional constraint on Sr 3 Ru207. 

When this work was complete, we have learned about a recent experimental observation of 
a large magnetoresistive anisotropy inside the dome-shaped transition line in S^R^Oy, 28] 
consistent with the existence of Fermi surface symmetry breaking in this compound. 

We are grateful to C. Honerkamp, A. P. Mackenzie, and R. S. Perry for valuable discus- 
sions, and especially to W. Metzner for a critical reading of the manuscript. 
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FIG. 1: The <iFSD phase diagram in the plane of magnetic field h and temperature T for fi = 1.0 
and g = 1.0; the transition is second order for high T (T c 2nd ) and first order for low T (T c lst ); end 
points of the second order line are tricritical points (T* ri ); the dotted line (h = 0.4) represents the 
van Hove energy of the up-spin band; the dashed line (T*) denotes a peak position of the uniform 
magnetic susceptibility and the thin dashed line represents T* in the absence of the dFSD. 




FIG. 2: The mean-field solution at low T for fi = 1.0 and g = 1.0. (a) h dependence of the order 
parameter; note that j] = + i]K (b) A metamagnetic transition due to the first order dFSD 
transition; the result for g = is also shown by a dotted line. FSs for g = 1 (solid line) and (gray 
line) at h = 0.3 (c) and 0.5 (d); the deformation of the inner FS in (d) is hardly visible. 
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FIG. 3: T dependence of \ f° r several choices of h for g = (a) and g = 1.0 (b); the dashed lines 
in (b) are data for g = 0; h = 0.4 corresponds to the van Hove energy of the up-spin band. 




FIG. 4: T dependence of 7 for several choices of h for g = (a) and g = 1.0 (b); a logarithmic T 
scale is used in (a); the dashed lines in (b) are data for g = 0. 
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